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Abstract. We consider a diffusion process witli coefficients tfiat arc periodic 
outside of an 'interface region' of finite thickness. The question investigated in 
this article is the limiting long time / large scale behaviour of such a process 
under diffusive rescaling. It is clear that outside of the interface, the limit- 
ing process must behave like Brownian motion, with diffusion matrices given 
by the standard theory of homogenization. The interesting behaviour there- 
fore occurs on the interface. Our main result is that the limiting process is a 
semimartingale whose bounded variation part is proportional to the local time 
spent on the interface. The proportionality vector can have non-zero compo- 
nents parallel to the interface, so that the limiting diffusion is not necessarily 
reversible. We also exhibit an explicit way of identifying its parameters in 
terms of the coefhcients of the original diffusion. 

Similarly to the one-dimensional case, our method of proof relies on the 
framework provided by Freidlin and Wentzell [FW93] for diffusion processes 
on a graph in order to identify the generator of the limiting process. 



1. Introduction 

The theory of periodic homogenization is by now extremely weU understood, 
see for example the monographs [BLP78, PS08]. Recall that the most basic result 
states that if X is a diffusion with smooth periodic coefficients, then the diffusively 
rescaled process X^{t) = eX{t/e^) converges in law to a Brownian motion with an 
explicitly computable diffusion matrix. If one considers diffusions that are 'locally 
periodic', but with slow modulations over spatial scales of order e~^, then it was 
shown in [BMP05] that the rescaled process converges in general to some diffusion 
process with a computable expression for both its drift and diffusion coefRcients. 

In this article, we will also consider the 'locally periodic' situation, but instead of 
considering slow modulations of the coefficients, we consider the case of a sharp (i.e. 
of size 0(1)) transition between two periodic structures. In the (much simpler) one- 
dimensional case, this model was previously studied in [HMIO], where we showed 
that the rescaled process converges in law to skew Brownian motion with an explicit 
expression for the skewness parameter. In higher dimensions, this model has not 
yet been studied to the best of our knowledge. The aim of this article is to clarify 
what is the behaviour of X'^ near the interface for very small values of e. It is 
important to remark at this stage that we do not make the assumption that our 
diffusion is reversible. As we will see in Section 2 below, there are then situations in 
which the limiting process is not reversible either, contrary to the one-dimensional 
situation. 

One feature of the problem at hand is that there is no finite invariant measure 
built into the framework of the problem. This is unlike most other homogenization 
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problems, even those exhibiting rather 'bad' ergodic properties, such as the ran- 
dom environment case [PV81, 01194] or the quenched convergence results for the 
Bouchaud trap model [BAC07]. Since in our case the invariant measure /i of X is 
only CT-finite, this leads to two problems when trying to compiite the effect of the 
behaviour of X near the interface in the limit e -4 0. Indcicd. one; would 'naively' 
expect that an effective drift along the interface can be described by the quantity 

y b{x)tJ,{dx) . (1.1) 

One problem with this expression is that there is no obvious natural normalisation 
for /i. Furthermore, since b is periodic away from the interface and the same is 
(approximately) true for yu, this integral certainly does not converge, even if we 
consider it as an integral over K x T"^"^ by making use of the periodic structure 
in the directions parallel to the interface. See however (2.4) and Proposition 6.3 
below for the correct way of interpreting (1.1) and our main result. Theorem 2.4 
below, on how this quantity appears in the construction of the limiting process. 

Another common feature of many homogenization results is the usage of a glob- 
ally defined corrector function to compensate for the singular terms appearing in the 
problem. This is of course the case for standard periodic homogenization [BLP78] , 
but also for a mimbcr of stochastic homogenization problems, as for example in 
[OS04, PV81, 01194, RhoOQ]. For the present problem however, it will be conve- 
nient to make use of corrector function that only cancels the singular terms away 
from the interface and to treat the behaviour of the limiting process at the interface 
by completely different means. 

One very recent homogenization result where discontinuous coefficients appear in 
the limiting equation can be found in [BEP09] (which in turn generalises [KKOl]). 
However, their framework is quite different to the one considered here and doesn't 
seem to encompass our problem. Much more closely related problems are homog- 
enization problems with the presence of a boundary [AA99, GVM08]. Those have 
been mostly studied by analytical tools so far. In our probabilistic language, what 
comes closest to the boundary layers studied in these articles is the cr-finite invariant 
measure of X , which is shown in Proposition 5.5 below to converge exponentially 
fast to a measure with periodic densities away from the interface. 

For simplicity, we will consider the case of a constant diffusion matrix, but it is 
straightforward to adapt the proofs to cover the case of non-constant diffusivity as 
well. More precisely, we consider the family of processes X^ taking values in R'^, 
solutions to the stochastic differential equations 

dX^ = ^b(J^^ds + dB{s) , X%0) = x, (1.2) 

where i? is a d-dimensional standard Wiener process. The drift b is assumed to be 
smooth and such that b{x + e^) = b{x) for the unit vectors Cj with i = 2, . . . ,d (but 
not for i = 1). Furthermore, we assume that there exist smooth vector fields b± 
with unit period in every direction and ?7 > such that 

b{x) = b+{x) , xi > T] , b{x) = b-{x) , xi < —rj . 

Figure 1 is a typical illustration of the type of vector fields that we have in mind. 

If we denote by X the same process, but with £ = 1, then the process X^ given 
by (1.2) is equal in law to the diffusive rescaling of X by a factor ^. In the sequel, 
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Figure 1 . Example of a vector field h satisfying our conditions. 



we denote the generator of X by £ and the generator of by We furthermore 
denote by L± the generators for the diffusion processes on the torus given by 

dX^ = h±{X^)ds + dB{s) , (1.3) 

and by /U± the corresponding invariant probability measures. With this notation 
at hand, we impose the centering condition b±{x) iJi±{x) = 0. 

Under these conditions, our main result formulated in Theorem 2.4 below states 
that the family X^ converges in law to a limiting process X. Furthermore, we 
give an explicit characterisation of X, both as the unique solution of a martingale 
problem with some explicitly given generator and as the solution of a stochastic 
differential equation involving a local time term on the interface {xi =0}. In 
addition to the homogenized diffusion coefficients on either side of the interface, 
this limiting process is characterised by a 'transmissivity coefffcient', as well as by 
a 'drift vector' pointing along the interface. 

The remainder of this article is structured as follows. After formulating our main 
results in Section 2, we show tightness of the family in Section 3. In Section 4, we 
then formulate the main tool used in the identification of the limiting process, 
namely a multidimensional analogue of the tool used by Freidlin and Wentzell in 
[FW93] to study homogenisation problems where the limiting process takes values 
in a graph. Section 5 is then devoted to the computation of the transmissivity 
coefficient, whereas Section 6 contains the computation of the drift vector. Finally, 
we show in Section 7 that the martingale problem is well-posed and we identify its 
solution with the solution to a stochastic differential equation. 

1.1. Notation. We define the 'interface' of width K by 

J'k = {x (^W^ : xi & [-K,K]} . 

We also denote by d^K its boundary. 

Frequently throughout the paper we will construct successive escape and sub- 
sequent reentry times particularly when constructing invariant measures in terms 
of the invariant measure of an embedded Markov chain as in [Has60]. We will de- 
note such pairs of stopping times as a, (j), which denote escape and reentry times 
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respectively. Other stopping tiraes not part of such a sequence will be denoted by 
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2. The Main Result 

Before stating the main result we will first define the various quantities involved 
and their relevance. It is clear that, in view of standard results from periodic 
homogenization [BLP78, PS08], any limiting process for should behave like 
Brownian motion on cither side of the interface = {xi = 0}, with effective 
diffusion tensors given by 

^ti = I i^ik + dkgt){Skj + dkgf) dn± . 

(Summation of k is implied.) Here, the corrector functions g±:f'^^ K** are the 
unique solutions to C±g± = —b± such that 



/ g±{x) li±{dx) = . 



Since h± are centered with respect to such functions do indeed exist. 

This justifies the introduction of a differential operator C on defined in two 
parts by £+ on 7+ = {xi > 0} and C- on 7_ = {xi < 0} with 



-^d.dj , (2.1) 



then one would expect any limiting process to solve a martingale problem associated 
to jC. However, the above definition of jC is not complete, since we did not specify 
any boundary condition at the interface J^q. 

One of the main ingredients in the analysis of the behaviour of the limiting 
process at the interface is the invariant measure fj, for the (original, not rescaled) 
process X. It is not clear a priori that such an invariant measure exists, since X is 
not expected to be recurrent in general. However, if we identify points that differ 
by integer multiples of ej for j = 2, . . . , d, we can interpret X as a process with 
state space R x T'^^^. It then follows from the results in [Has60] that this process 
admits a cr-finite invariant measure on M x T''"-'^. 

Note that the invariant measure /x is not finite and can therefore not be nor- 
malised in a canonical way. However, if we define the 'unit cells' by 

C; = [j, j + l]x T'^-i , Cr = [-j - 1, -j] X T-^-i , 

then it is possible to make sense of the quantity q± = Hmj^^ fJ-i^f) (we will show 
in Proposition 5.5 below that this limit actually exists). 
Let now p± be given by 

^ q±Dti 

q+D++q-D^,' 
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which can also we rewritten in a more suggestive way as 

P+ _ q+Dti 



(2.2) 



p- q-D-^^i 

This is the homogenized diffusion coeffieient in the direction perpendicular to the 
interface, weighted by the invariant measure of a unit cell. Comparing with the 
one- dimensional case [HMIO], one would expect this to yield the likelihood for 
to exit a small (but still much larger than e) neighborhood of the interface on a 
specific side. 

Remark 2.1. The ratio 

(2.3) 



gives the asymptotic probability of the process being located in the rhs (+) of the 
interface after a long time. This follows from the weak convergence of the first 
component to a skew Brownian motion with (possibly) different diffusion coefficients 
on either side of the interface. If we rescale this skew BM on either side of the 



interface by \ Dfi to obtain a standard skew BM, we can use the scale function of 

BM to finish the verification of (2.3). 

However, unlike in the one-dimensional case, these quantities are not yet suffi- 
cient to characterise the limiting process. The reason is that since is expected 
to spend time proportional to e in the interface, but the drift is of order there, 
it is not impossible that the limiting process picks up a non-trivial drift along the 
interface. It turns out that this drift can be described by the coefficients aj given 
by 

= + ^) / (^^- W + ^9j{x)) '"(^^) ' (2-4) 

where n is again normalised in such a way that + q- = 1 and where g is any 
smooth function agreeing with g± on either side of the interface (see Section 3 
below). 

Remark 2.2. Since /j^^-jd-i C(p{x) ii{dx) = for every smooth compactly supported 
function (j), one should interpret the integral on the right hand side of (2.4) as a 
'renormalised' form of the intuitive more meaningful quantity (1.1). 

Remark 2.3. The expression (2.4) is useful in order to generate examples with 

non-vanishing values for the coefficients Ui. 

Given all of these ingredients, we can construct an operator C as follows. The 
domain S>{C.) of C consists of functions / : M** — )• R such that 

• The restrictions of / to /+, and are smooth. 

• The partial derivatives dif are continuous for i > 2. 

• The partial derivative dif{x) has right and left limits d\f\i^ a,s x ^ J{j 
and these limits satisfy the gluing condition 

d 

P+difW - P-dif\i- + "^odjf = . (2.5) 

j=2 
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Figure 2. Sample paths at small (left) and large (right) scales. 

For any / e S>{C), we then set Cf{x) = C±f{x) for x G I±. With these definitions 
at hand, we can state the main result of the article: 

Theorem 2.4. The family of processes converges in law to the unique solution 
X to the martingale problem given by the operator L. Furthermore, there exist 
matrices M± and a vector K gR''' such that this solution solves the SDE 

dX{t) = lx^^oM_dW{t) + lx,>oM+dW{t) + K dL{t) . (2.6) 

where L denotes the symmetric local time of Xi at the origin and W is a standard 
d-dimensional Wiener process. The matrices M± and the vector K satisfy 

M±Ml = D^ , Ki=p+-p_ , Kj = aj , 

forj = {2,...,d}. 

In Figure 2, we show an example of a numerical simulation of the process studied 
in this article. The figure on the left shows the small-scale structure (the periodic 
structure of the drift is drawn as a grid). One can clearly see the periodic structure 
of the sample path, especially to the left of the interface. One can also see that 
the effective diffusivity is not necessarily proportional to the identity. In this case, 
to the left of the interface, the process diffuses much more easily horizontally than 
vertically. 

The picture to the right shows a simulation of the process at a much larger scale. 
We used a slightly different vector field for the drift in order to obtain a simulation 
that shows clearly the strong drift experienced by the process when it hits the 
interface. 

Remark 2.5. Since the quadratic variation of X has a discontinuity at Xi = 0, 
we do have to specify which kind of local time L is. Using the symmetric local 
time yields nicer expressions. See for example [RY91, Lej06] for a definition of the 

symmetric local time. 

Analyzing what this means for a simple example, we consider the case of a two 
dimensional problem where we have 6i = and 62 = f{xi) for / a smooth function 
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that is zero outside of Jf,. Clearly p± = ^. In this case the invariant measure /i of 
the process X is given by 5 times Lebesgue measure on R x 5^ and we can choose 
g = 0. This implies that we then simply have 




as one would expect. 

3. Tightness of the family 
The aim of this section is to prove the following tightness result: 
Theorem 3.1. Denote by W the law of on C(M+,R''). Then the family 

Similar to what happens in the classical theory of periodic homogcnization, it 
will be very convenient to construct a 'corrected process' Y, obtained by adding 
to X a corrector function that cancels out to first order the effect of the small 
oscillations. To this aim, wc introduce a smooth function gilS.'^ — > M'' which is 
periodic in the directions 2,...,c? and such that g{x) = g+{x) for xi > rj and 
similarly for Xi < (Recall that g± were defined in Section 2.) We do not 

specify the behaviour of g inside the interface except that it has to be smooth 
in the whole space and periodic in the directions parallel to the interface. We fix 
such a function g once and for all from now on. Wc furthermore denote by the 
process defined by = X^ + eg{e'^X^), as well a.^ y — x + eg{x/e) for its initial 
condition. 

Defining the corrected drift b{x) = (Cg + 6) (x) and the corrected diffusion coef- 
ficient d'ij{x) = 6ij + djgi{x), it follows from Ito's formula that the zth component 
of Yy satisfies 

(r;) .(t) = yi + h.{^^xi{s)) ds + £ {}^K{s)) dw^{s) . (3.1) 

It is very important to note that the corrected drift h vanishes outside of J^,, so 
that the process Y is subject to a large drift only when X is inside the interface. 

Our main tool in the proof of Theorem 3.1 is the following result, which is very 
similar to [SV79, Thm 1.4.6]: 

Proposition 3.2. Let !^ he a family of probability measures on Q. = C(]R+,M'') 
and denote by x the canonical process on fi. Assume that 

lim sup P(|a;(0)| > i?) = . 

Furtherm,ore, for any given p > 0, let tq = 0, and define recursively Tj+i = 
infj>^^ \x{t) — x(Ti)\ > p. Assume that the limit 

lim esssupP [t„+i — Tn < <5|>^t„] —^0, P a.s., on {t„ < ooj (3.2) 

(5->0 

holds uniformly for every P e ^ and every n > 0. Then the family of probability 

measures ^ is tight on f2. 

Proof. The proof is similar to that of Theorem 1.4.6 in [SV79], except that their 
Lemma 1.4.4 is replaced by (3.2). 
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Fix an arbitrary final time T > 0. Furthermore, denote for w € f2, 

Np = Np{ui) = min{n : r„+i > T} . 

and tlie modulus of continuity by 5p, 

Sp = Sp{uj) = min{r„ - r„_i : 1 < n < Np{w)} . 

Note that this expression depends on p via the definition of the stopping times r,. 

With this notation at hand, tightness follows as in [SV79] if one can show that 
lim5_>o suppg^ P((5p < (5) = for every fixed p > 0. As in [SV79], one has for every 
fc > the bound 

k 

< <5) < ^E[P -n<5\^rA] + n^p > k) . 

i=l 

For every fixed fc > 0, the first term then converges uniformly to by assumption. 
Since the second term is independent of 5, it remains to verify that converges to 
as A; — >■ 00, uniformly over (convergence for every fixed P e is trivial but not 
sufficient for our needs). 

This is a consequence of [SV79, Lemma 1.4.5], provided that one can find A < 1 
such that E[e-('"'+i-^')| < A. This in turn follows from 

E[e-(ri+i-T,)| < P [t,+i -n< to\^r^] + e-*°P [n+i -n> to\ 

< e-*° + (1 - e-*«)P [Ti+i - Ti < toi =^n] • 

Indeed, by choosing sufficiently small, this term can be made strictly less than 

1, provided that P [t^+i — < to| -^tJ tends to zero uniformly (over the members 
of ^ and over i) as to tends to zero, which is precisely our assumption. □ 

We now turn to the 

Proof of Theorem 3.1. Recall that we defined the process = X'^ + eg{e~^ X'^) in 
Section 2. Note then that, just as in [HMIO, Prop. 2.5], the tightness of the laws of 
X^ is equivalent to that of the laws of Y^. Therefore, all that remains to be shown 
is that we have the bound (3.2) for the law of , uniformly over e G (0, 1]. The 
approach that we use is to consider separately the martingale part and the bounded 
variation part for Yy given by (3.1), and to show that the probability of either of 
these moving by at least | during a time interval 5 tends to zero uniformly over 
the initial condition. 

Given any fixed p,7 > 0, we want to show that there exists a sufficiently small 
5 > Q such that P(r„+i - t„ < 5\^r„) < 7 uniformly over P G ^ (that is 
uniformly over the laws of Y^ with e G (0, 1]) and n. We split the contributions 
from the martingale and the bounded variation parts in the following way: 



P(t„+1 -Tn<S\^rJ= Px(r„) ( SUp \Y{t) - Y{0)\ > p) 

^ t<S ' 

<supPa.(sup i / hi{e-^Xl{s))ds > J) 
+ supP^(sup \ u^{£-^Xl{s))dWj{s 

X ^ t<<5 Jo 



2 



< — supE^ / |6i(e-iX|(s))|rfs 

X Jo 
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+ - sup Ea; sup 
P X t<5 



faij{e-'X^{8))dWj{s) . (3.3) 
Jo 



Here, we used the Chebychev's inequality to obtain the last bound. Since the 
functions aij are uniformly bounded, the stochastic integral appearing in the second 
term is easily bounded by 0{VS) by the Burkholder-Davis-Gundy inequalities. 
Furthermore, by the definition of the corrector function g, there exists ?) > such 
that b{x) = for a; J^^, so that there exists a constant C such that 

P(T„+i-r„ <5|^,J < — supeJ/ i^.^(x|(s))ds) +^ . (3.4) 

For fixed p > 0, the second term obviously goes to as (5 —> 0, uniformly in £, so it 
remains to consider the first term. As one would expect from the expression for the 
local time of a Brownian motion, it turns out that the expected time spent by the 
process in J^^ scales like e\fb, thus showing that this term is also of order ^/S/ p. 
Once we are able to show this, the proof is complete. 

The occupation time of the interface appearing in the first term of (3.4) is 
bounded by the trivial estimate Cd/{pe), which goes to as (5 — )■ provided that 
we consider e > VS, say. We can therefore assume without any loss of generality 
in the sequel that wc consider s < VS. 

The idea to bound the occupation time is the following. We decompose the 
trajectory for the process into excursions away from the interface, separated 
by pieces of trajectory inside the interface. We first show that if the process starts 
inside the interface, then the expected time spent in the interface before making a 
new excursion is of order e^. Then, we show that each excursion has a probability 
at least e/VS of being of length 5 or more. This shows that in the time interval 
S of interest, the process will perform at most of the order of VS/e excursions, so 
that the total time spent in the interface is indeed of the order eVS, thus showing 
that the first term in (3.4) bchavcis like VS/p, as expected. 

More precisely, we first choose two constants K > and K > such that the 
chain of implications 

{X' e 4e} ^ {Y' e J ^ {X' e ^k-ds} ^ {X' e j^ks} , (3.5) 

holds. We then set up a sequence of stopping times in the following way. We set 
(f)o = and we set recursively 

CTn = inf{t>0„ : X^t) (^.fKe) . 

<t>n = inf{t > On-i ■■ Y'{t) e ^^J . 

(Note that we can have (Tq = if the initial condition does not belong to J^£. Apart 
from that, the second implication in (3.5) shows that increments from one stopping 
time to the next are always strictly positive.) This construction was chosen in such 
a way that the times when X^ € J^e always fall between ^„ and cr„ for some n > 0. 
In particular, if we set 

N = inf{n > : (^„+i - cr„ > (5} , 

then we have the bound 

1^,. (X'As)) ds) < supEj^(a„ - 0„)) 

^ n=0 
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= SUp^Ea;((cr„ - ^„)ljv>n) = ^SUpPa;(A'' > n) SUp E^; (Exs(0„)O-i) 

^ n=0 n=0 ^ ^ 

where we used the strong Markov property and the fact that {N > n} is 
nicasurablc in order to obtain the last identity. It follows from the definition of A'' 
that this expression is in turn bounded by 

sup E^cTo > ( sup P^(0o < S)) = — 7 " . 

We now bound both terms appearing in this expression separately. 

First, we turn to the expected escape time from the interface, ExO'o- The idea 
is to use a comparison argument just like in [HMIO, Proposition 3.8]. We define 
a 'worst-case scenario' process V^, which is the solution to the SDE with initial 
condition x, diffusion coefficient 1 and drift coefficient given by by, where 



for some constant by > 0. We then have: 



for a; > 0, 
for a; < 0, 



Lemma 3.3. There exist by > and K > such that, if we define = mf{t > 
: V^{t) ^ J^Ke}' "'e have 

E^ao < E^r* , 

for every a; e M''. 

The proof of Lemma 3.3 is almost identical to that of [HMIO, Proposition 3.8], 

so we are going to omit it. A straightforward calculation using the particular form 
of the drift coefficient for V allows to check that there exists indeed a constant 
C > such that the bound 

supE^r^ < Ce^ , 

X 

holds so that, combining this with Lemma 3.3, we have sup^g^d ExCTq < Ce^ . 

Let us now turn to the bound on fx{4'o > The idea here is to look at 
the process instead of and to time-change it in such a way that we can 
compare it to a standard Brownian motion. Note first that the last two implications 
in (3.5) show that if we start with anywhere outside of J'ks, then the first 
component of has to travel by at least e before the process Y^ can hit J^g- 
Furthermore, it follows from (3.1) that the time change Ct such that Y'^{Ct) is a 
standard Brownian motion satisfies Ct > ct for some c > 0. It therefore follows 
that, setting H{z) = inf(>o{i?( > z}, one has the lower bound 

inf P^(0o >5)> V{H{e) > 6/c) . 

The explicit expression for the law of B.[z) given in [BS96, p. 163, eq. 2.02] yields 
in turn 

,oo ^-l/(2t) 

P(i/(£) > 6/c) = / 3 dt . 



It follows immediately that this in turn is bounded from below by Ce/\/^ for some 
C > 0, provided that e < VS. Collecting these bounds completes the proof of 
Theorem 3.1. □ 
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4. Main tool for identifying the limit process 



Instead of considering a graph as before, we will consider a generalized multidi- 
mensional version different from that considered by Freidlin and Wentzell in [FW93, 
Section 6]. Note that the generalisation considered here is different (and actually 
simpler) than the one considered in [FW06] . We consider processes in R'^ and we set 
7_ = {a; € M'' : xi < 0}, and similarly for /+. We consider a family of R''-valued 
processes X'^ and we denote by the first hitting time of Correspondingly, 

is the first escape time of the set by X'^. 

With this the main tool will be the following multidimensional analogue of 
[FW93, Thm 4.1]: 

Theorem 4.1. Let Ci he second order differential operators on li with bounded 
coefficients and let Di be some sets of test functions over li whose members are 
bounded and have bounded derivatives of all orders. Suppose that for i € {+,—}; 
any function f £ Di and for any A > 0, the bound 



e-^^'f{X%T'))-f{X%0)) 

+ f e-^'(\f{X\t))~Cif{X'{t)))dt 



0{k{e)) 



(4.1) 



holds as e — > 0, uniformly with respect to x G li. Assume furthermore that the rate 
k is such that lim^^o k{e) = 0. 

Assume thai, for every A > and every i G {+, — }, there exist functions Ui^x € 
Di such that CiUi^\{x) = Xui^\{x) holds for x & li with \xi\ < 1 and such that 
u±^\{x) — 1 for xi = and xi = ±1. 

Assume that there exists a rate 5 = 5{e) — )• such that 5{e)/k{e) oo as £ ^ 
and such that for A > 0, 



E! 



.S,5){X!{t)) dt 







(4.2) 



as e ^ 0, uniformly in the initial point. Assume the convergence 

F%[X^{T')Gli]^Pi, (4.3) 

holds uniformly in x in the set J^^ for some constants p± withp^+p- = 1. Assume 
furthermore that there exist constants aj and C such that 



X|(r^) 



a 



— E^ 
6^ ' 



{X^{r')-x,y 



<C 



(4.4) 



for j > 2. Again, the limit is assumed to be uniform over x G .^r; as e — ^ 0, and 
the inequality is assumed to be uniform over all e G (0, 1] and all x G J^,,. 

Let then D be the set of continuous functions / : R"^ — >■ R such that the restriction 
of f to li belongs to Di and such that the gluing condition (2.5) holds. Then, for 
any fixed f & D, to >0 and A > 0, 



A(£) = ess sup 



ES 



-At 



to 



Xf{X%t))-£f{X%t)) 



dt 



— e 



-At, 



°f(.X^to)) 



[0,to] 







(4.5) 



as £ — >■ 0, uniformly with respect to x. In particular, every weak limit of X'^ as 
£ ^ satisfies the martingale problem for C. 
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Remark 4.2. Note that we did not specify how 'large' the sets Di of admissible test 
functions need to be. If these sets are too small then the theorem still holds, but 
the corresponding martingale problem might become ill-posed. 

Proof. Since the proof is virtually identical to that of [FW93, Thm 4.1], we only 
sketch it here. The basic idea behind the proof given by Freidlin and Wentzell 
is to rewrite (4.5) using the Strong Markov Property of as a sum of terms 
between successive stopping times. To this effect, set for example ctq = and then 
recursively = inf{t > ct„ : Xl{t) e J^en}, cr„+i = \ni{t > (pn ■ Xf(t) ^ Ji}. 
They then break up the term produced from (4.5) into two sums of analogous terms 
between times an and (j)n and those between (f>n and cr„+i. 

The terms covering the time intervals [cr„, 0„] are bounded exactly as in [FW93], 
making use of (4.1), together with the bound X)™ -'^^^e"'^'^" = 0{1/S) which follows 
from the existence of the functions Ui^x just as in [FW93]. 

Using assumption 4.2, the terms covering the time intervcils [^^i , c7„+i] are then 
simplified to 

n 

modulo contributions that converge to as £ 0. Since the expectation of this 
term is bounded by 

sup E4f{X%T')) - /(x)) ^Ee-^<^" , 

and since we already know that X^n^^"'^*" ^ ^(V*^); remains to show that the 
supremum is of order o{S). It follows from Taylor's expansion and the fact that 

f GC^ outside of the interface, that on the event f2+ =' {Xf (r"^) > 0}, one has 

d 

f{X'{r')) - fix) = 6dif{x)\i^ +J29if{x){Xt{r') - x^) + (r*) - Xi\') , 

i=2 

and similarly on il^ = {Xf{T^) < 0}. Combining this with (4.4), we thus have 
E,{f{X'{T^)) - f{x)) = 6dif{x)\i^^,{^+) + 8dif{x)\i_V,{n-) 

d 

+ 5Y,Oiidif{x) + o{5) . 

Since we assume that ¥x{^±) — >■ p± uniformly over x S J^e, the required bound 
now follows from the gluing condition. □ 

Most of the remainder of this article is devoted to the verification of the assump- 
tions of Theorem 4.1. The bounds (4.1) and (4.2) will be relatively straightforward 
to verify and this will form the content of the remainder of this section. The con- 
vergence (4.3) is the one that is most difficult to obtain and will be the content 
of Section 5. Finally, we will show that (4.4) holds in Section 6. We start by the 
following result: 

Lemma 4.3. Let C± be as in (2.1) and let X^ be the family of processes from 
Section 2. Then, the bound (4-1) holds with k{e) = e for every A > and for every 
smooth bounded function f-.Ii^R that has bounded derivatives of all orders. 
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Proof. It follows from [HMIO, Lemma 3.4] that, for any initial point x with xi 
and for e sufficiently small so that x J^,,, 



0(e) , (4.6) 



for h centered with respect to /x+ (resp. fi- if xi < 0). We assume that a;i > 
from now on, but the calculations are identical for the case xi < 0. 

Note now that it suffices to obtain the bound (4.1) for the family of processes , 
since \\Y^{t)-X^{t)\\ = 0(e), uniformly. Applying Ito's formula to e-^'"V(^^(T"^)), 
we obtain the identity 

e-^^'f{Y%Tn)=f{y)+ r -Xe-^^f{Y%s))ds 

Jo 

e-'''a,J^^)d^f{Y%s))dWk{.s) . 



Since — X^l < 0{e) and since all derivatives of / are assumed to be bounded, 
it then follows from (4.6) that 

E(e-^-V(y^(r-))) =/(?/)- AE J' e-^' f{Y%s)) ds 

1 

+ -E e-^^D+dlf{Y-{s)) ds + 0{e) , 

which is precisely the required result. □ 

Additionally wc have that the solution to CiU = Am on 7^, u = 1 on {xi = 0} and 
= ±1}, is bounded and has bounded derivatives of all orders. This follows from 

the fact that u is given explicitly by u{x) = Cie'v^'^^^"^"^^! + C2e~ v'^^^^^"^^! for 
some constants Ci. We now show that the process Y^ satisfies the bound (4.2), i.e. 
it does not spend too much time in the vicinity of the interface: 

Lemma 4.4. // we choose 5 ~ e" for any a & (|, 1), then (4.-2) holds for the family 

of processes from Section 2. 

Proof. Again, it suffices to show the bound for the process Y^ since it differs from 
X^ by 0{e). We would like to use an argument similar to what can be used in the 
one-dimensional case [HMIO], that is we time-change the corrected process Y"^ in 
such a way that it becomes a diffusion with diffusion coefficient 1. Its drift then 
vanishes outside of the interface and is bounded by K/e for some K > Q. At this 
stage, one compares this process to the 'worst-case scenario' process given by 

dZ^ = b{Z') dt + dB{t) , 

where the drift b is given by 

-/sTe-i ifze[0,Zc), 
b{z) = { Ke-^ ifze(-/£,0), 
otherwise. 
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for some / e R. It can then be shown that spends more time in the interface than 
does, so that the requested bound can be obtained from a simple calculation. 
The problem with this argument is that in the multi-dimensional case the time- 
change required to turn the first component of into a diffusion with unit diffusion 
coefficient is given by 

Tt = infjs e M+ : / X^((5ii + d,gi{e-^X'{u))^ du > (4.7) 

"'O i=l 

We do not know of an argument giving a uniform bound from below on the quantity 
appearing under the integral in this expression. Therefore, an upper bound on the 
time spent by the process in the interval {—S, S) does not give us any control on 
the time spent by Y'^ (and therefore X^) in that interval. 

Because of this, we modify our argument in the following way. We break up the 
integral in (4.2) as 



(-5,5) 



{yiit))dt 



= E, 









( — c£, ce) 



(Yat)) dt 



+ E, 



[f 

Jo 

f 

Jo 



-\t-l 



.S,-oe){Ynt)) dt 



(4.8) 
(4.9) 



-Af 1 



L(c£,5) 



{Ynt))dt 



where Yf is the first component of Y^ and c is a value to be determined. By 
symmetry, the last two terms are of the same order, so that it is sufficient to bound 
the first two terms. In order to bound the first term, we use the argument outlined 

above, but we replace F*^ by the process Y^ given by Y^{t) = X^ (t) + sg{e~^ X^ (t)) , 
where the corrector g has the following properties: 

1. The function g{x) is smooth, periodic in the variables parallel to the inter- 
face, and equal to g{x) for x ^ J^^ for some ci. 

2. One has the implication Y^ G J^e => Y^ G J^^^ for some C2 < Ci. 

3. If Y'' e ^^26, then g{e-^X^) = 0. 

It is always possible to satisfy these properties by choosing ci sufficiently large 
and setting g = in a sufficiently wide band around the interface. We now set 
Z{t) = Y{Tt), where Tt is defined as in (4.7), but with g replaced by g, so that it 
follows from the second property that one has the bound 



f° 

^-cs,cs){Yat))dt <E, / 
Jo 



<E, 



i-C2e,C2e){Ynt))dt 

(-C2£,C2£)(2'l(i)) dTt . 



At this stage, we remark that since the function g has bounded derivatives, there 
exists a constant Ki such that Tt > Kit almost surely. On the other hand, it 
follows from the last property that one actually has dTt = dt whenever Y^ € ^C2£i 
so that this expression is bounded by 



E, 



/>CXD 

.1 e-^^*l(_,,,,,,,)(Zi^(i))dt. 
Jo 

This expression in turn can be bounded by 0{e) just as in [HMIO]. 
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We now proceed to bounding the term (4.9). For this, let us first introduce a 
constant C3 < c and make c from (4.8) sufficiently large such that 

4. The impHcation X^t) G J^cse Y^{t) e J^cs holds. 

5. One has C3 > rj+1. 

Then, we define a series of stopping times {(l)'n}n and {cr^}n recursively by (p'_i = 
0, . . ., ff,; = M{t > <_i : Xf{t) {-25, -cse+e)} and < = inf{i > a'^ : Xf{t) e 

(-5,-C3£)}. 

Now we can use the Strong Markov property as in [FW93, Lemma 4.1] with 
the stopping times to obtain the bound Ej, [^^q e^'^'^"^^^] = O(-), uniformly 
in the initial point x for x € {x : Xi = —Cse + s} U {x : xi = —26}. This is 
a consequence of the fact that Ea;[e~'^'^o] = 1 — 0(e) uniformly. Furthermore, it 
follows from the definition of these stopping times, property 4., and the strong 
Markov property that (4.9) is bounded by 

/ e-^*l(_5,_e3e) (xm) dt<E,J2 / e-^* dt 

< X-'E^ Yl e"^^"-^ K - 'P'n-i) (4-10) 

n>0 

00 „ 

< A-i(e^ Ve-^<(^)) supE^CT^, < — supE^a^, . 

^ ^ ^ ^ ^-^ ^ 

It follows that it suffices to be able to choose 6 in such a way that E^ctq is o(e) 
uniformly in the initial point. Specifically, we will show that (4.10) is 0{6^), so 
that the claim follows. 

This will be a consequence of the following result: 

Lemma 4.5. Let X~ be as in (1.3) and define X~'^{t) = sX~ {e~^t). Let r = 
'mi{t > : X^'^{t) ^ [—1,0]}. Then, there exists a constant C such that 

ExT < C , 

independently of e G (0, 1] and independently of x gR"^. 

Before we prove Lemma 4.5, we use it to complete the proof of Lemma 4.4. It 
follows from property 5. that up to time ctq, the process X^ is identical in law to 
the process Furthermore, the stopping time ctq is certainly bounded from 

above by the first exit time of the first component of X^-^ from {—26, 0). Rescaling 
space by a factor 26 and rescaling time correspondingly by 46^, we deduce from 
Lemma 4.5 that Ectq < AC 6^, uniformly in the initial condition as required. □ 

Wc now turn to the 

Proof of Lemma 4-5. Denote by U the region {x G R'^ : xi G [—1,0]} and define 
f^ by f^{x) = E^T. Then satisfies 

C'f' = -1, f%x) = hi X GdU, 

where = + e~^b-{e~^-)Vx- In order to obtain a bound on /, we will give a 
uniformly bounded (uniformly over e) function such that it satisfies, 

£V = -1, f{x) > for a; e dU. (4.11) 
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C^g^ix) = -Ci - (J) (1 + 2xi) . 



It then follows from the strong maximum principle (which we can apply since our 
diffusion is periodic in the directions in which U is unbounded) that > f^, so 
that the requested bound holds. 

We use a standard multiscale expansion for g^ of the form 

= .90 + £.91 + £^32 • 

Now to find such a g^. We proceed by starting off with a constant order term, that 
is, the typical term one would expect for the escape time if we were dealing with a 
Brownian motion, then removing the order ^ terms that arise when the operator 
acts on the constant order term by adding an order e term. Then finally, we 
add an order term to remove the constant order terms that are produced by the 
action of on the order e term. Incidentally, this approach of correction works 
exactly with the maximum order term in e being 2 and produces a series of terms 
that are known and have the right properties to provide a uniform bound. 

Taking guidance from the fact that the homogenized process is given by Brownian 
motion, we make the ansatz go{x) = C2 — CiXi{l+xi), for Ci and C2 two constants 
to be determined. Applying jC^ to go yields 

e 

for the first component of 6_. Our aim now is to choose gi in such a way that 
Cg\ contains a term of order that precisely cancels out the second term in this 
expression. Denote as in the introduction by g- the unique centered solution to 
the Poisson equation 

Cg. = 6_ , (4.12) 

where £ = + b-Vx is the generator for the non-rescaled process. We then set 
gi{x) = Ci(l + 2xi)g-^i{e~^x), where gf_^i is the first component of fif_, and we 
note that 

= — &-.if-)(l + 2,Ti) + Cii^(-) , (4.13) 

for some periodic function F independent of £ and of Ci. The term involving F 
appearing in this expression is still of order one, so we aim to compensate it by a 
judicious choice of 52- It is not necessarily centred with respect to the invariant 
measure of our process, but there exists a periodic centred function h such that 

Ch = F-K , K = j F{x) ii{dx) =- j \^g-^i{x)\^ n{dx) + 2 J fi{dx) . 

Finally, setting 52(2;) = —h{e^^x), wc obtain 

C'g' = Ci{K-l) = -Ci J |ei - Vg-,i{x)\^ fi{dx) . (4.14) 

Since the integral is strictly positive, the right hand side can be made to be equal 
to —1. Furthermore, since the corrector terms egi + £^(72 are uniformly bounded 
for e < 1, it is straightforward to find a constant C2 that ensures that g{x) > for 
X e dU, thus concluding the proof. □ 
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5. Computation of the transmissivity coefficient 
The aim of this section is to prove that 

Proposition 5.1. The identity (4-3) holds for the family of processes in Sec- 
tion 2 withp± given by (2.2). 

Let us first introduce some notation. Given a starting point x G J^,, we set 
_ f>^(^x{T'^^'>) > 0), and similarly for p!!'*^, where r*^*^' is the first hitting time 
of d^. We furthermore set 

-fc x,k k ■ c x.k (fe) I k\ 

p+ = sup p+' , pl = mf p+ , p\' ^ -{p++p1) , 

and similarly for p_. It is clear that Proposition 5.1 follows if we can show that p\ 
converges to a limit satisfying (2.2) and — p^ — )■ as — >■ cxd. 

We will first show the latter, as it is relatively straightforward to show. In 
order to show the convergence of p'ji , our main ingredient will be to show that the 
invariant measure ii{dx) for the process X looks more and more similar to ii±{dx) 
as xi — >■ ±00. Note that in this whole section, we will always consider X and X± as 
processes on M x T'^"-'^, obtained by identifying points (.t, y) such that xi = yi and 

~ Uj ^ ^ for j > 2. With this interpretation, the interface is compact and we 
will show that the processes are recurrent. If we were to consider them as processes 
in R'', they would not be recurrent for d > 3. 

Before we show that indeed p^ — p^ — )■ 0, we obtain some recurrence properties 
of X and ensure that it visits any open set in sufficiently often before the hitting 
time T^''^ . 

Lemma 5.2. Fix a neighborhood 7 C Then the probability for X to enter 7 
before hitting dJ^, starting from an arbitrary initial point in tends to 1 uniformly 
as k ^ 00. In particular, the process X is recurrent. 

Our first step in showing this result is to argue that if the process starts at 
distance 0(1) of the interface, then it will return to the interface with overwhelming 
probability before exiting J^: 

Lemma 5.3. There exists K > such that the probability, starting at x, for X to 
return to .f^ before hitting dJ)^, is bounded from above by 1 — and from below 
by 1^1^. 

Proof. Denote by /'^(a;) the probability of hitting before 9J^, starting from x. 
We assume without loss of generality that Xi > 0, since the case Xi < follows 
using the same argument. The function /'^ then satisfies the equation Cf'^ = 0, 
endowed with the boundary conditions f''{x) = 1 if Xi = and f'^{x) = if xi = fc. 
As in the proof of Lemma 4.5, we aim to construct a function g'^ satisfying Cg'' = 
and such that either g'^ix) < f^{x) on the two boundaries or g^{x) > f'^{x) on the 
two boundaries. The claim then follows from the maximum principle. 
Let g+ be as in (4.12) and set 

g''{x) = l-k-^{K + xi-g+,i{x)) , 

for some constant K to be determined. It is straightforward to check that g'^ 

docs indeed satisfy Cg'' = 0, as well as the required inequalities on the boundary, 
provided that K is either sufficiently large or sufficiently small. This concludes the 
proof. □ 
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We now use the result of lemma 5.3 to prove lemma 5.2. This is done using the 
strong Markov property in conjunction with success/failure trials. 

Proof of Lemma 5.2. Consider the two hypcrplancs that delimit J^q and two further 

hypcrplancs at distance m from with m a sufficiently large constant to be 
determined later. Wo then break the process into excursions from to d^ri+m 
and back. 

More precisely, we define two sets of stopping times {cr™},i and {ct>J" },i recursively 
by = mi{t> : X{t) G ... , C = inf{t > < ■ X{t) G ^J, 

(7^1 = inf{t > (f)^ : X{t) e dJ!r}+m}- We furthermore denote by the u-algebra 
generated by trajectories of X up to the time 0™ and by ^„ the cr-algebra generated 
by trajectories of X up to the time crJJ!|_i. We also denote by the first hitting 
time of the set 7 and by t^''^ the first hitting time of the set 9J^. 

It follows from the ellipticity of X and the resulting smoothness of its transition 
probabilities that there exists some p > such that infa^gaj^r^ ^1(2;, 7) = 2p > 0. 
Furthermore, it is straightforward, for instance using a comparison argument with 
a process with constant drift away from the interface and using the continuity of 
paths, to show that 

lim sup Px(aJ" < 1) = . (5.1) 

It follows that we can choose m large enough so that the probability appearing in 
(5.1) is bounded above by p. As a consequence, for such a choice of m, one has the 
almost sure bound 

P(r^ < C+i \^n)>P- (5.2) 
On the other hand, it follows from Lemma 5.3 that the probability that the process 
hits dJ^k between a™ and (f)^ is bounded from above uniformly by = 0{k~^) so 
that, almost surely, 

<C+il^n) </3fc. (5.3) 
Note furthermore that by construction the event appearing in (5.2) is ,^„-measurable. 

Denote now by ¥„ a Markov chain with states { — 1,0,1} such that {±1} are 
absorbing and such that P(r„+i -1 1 y„ = 0) = p, P(y„+i = 1\Y„ = 0) = Pk- 
As a consequence of (5.2) and (5.3), it is then possible to couple Y and X in such 
a way that the following two implications hold almost surely: 

{y„ = and y„+i = -1} {C < < <yn+i < r^"''} 

{<Vi < < C+1 < r^} => {{Yn - and r„+i = 1)} 

It follows that the probability of entering 7 before the hitting time t^''^ is bounded 
from below by 

P(r^ < T^^^h > P( lim y„ = -1) = . 

n->-oo P + Pk 

Since p is fixed and /3fe = 0{k~^), this quantity can be made arbitrarily close to 1. 

This shows that the set 7 is recurrent for X . Since furthermore X has transition 
probabilities that have strictly positive densities with respect to Lebesgue measure 
(as a consequence of the ellipticity of the equations describing it), recurrence follows 
from [MT93, Theorem 8.0.1]. □ 

We now use this result to prove 

Proposition 5.4. — — )• as k ^ 00. 
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Proof. The idea is to use; the faet that, before the proeess exits ,J^, it has had 
sufficient amount of time; to forget about its initial condition by visiting a small set 
on which a strong minorising condition holds for its transition probabilities. 
Fix a value /9 > 0. Our aim is to show that there then exists ko > such that 

say, for every k > ko. Since p^'^ = 1 — p^'', the claim then follows. We restrict 
ourselves to the bound for p+ since the other bound can be obtained in exactly the 
same way. 

The argument is now the following. It follows from the smoothness of transi- 
tion probabilities that there exists a neighbourhood 7 of the origin such that the 
transition probabilities at time 1 for X, starting from 7 satisfy the lower bound 

p{y) = int Pi{x,y) , 

with /jj^^d-i piy) dy > 1 — f3/2. It then follows immediately that for a; G 7, one has 
^x^k ^ ^(hfc _ -f^{3t<l : X{t) e dJi). For arbitrary x, it therefore follows 
from the strong Markov property that 

^x,k ^ pO,k _ _ supPj,(3i < 1 : X{t) &Jk)- Px(X hits dJu before 7) . 

j/67 

The last term can be made smaller than /3/4 by Lemma 5.2. The remaining term 
Pj,(3f < 1 : X{t) e on the other hand was already shown to be arbitrary small 
in (5.1). □ 

We next show that the invariant measure of the process converges to that of the 
relevant periodic process with increasing distance from the interface. 

Proposition 5.5. Let A denote a hounded measurable set and denote by fi the 
(unique up to scaling) invariant a-finite measure of the process X. Denote fur- 
thermore by fx± the invariant m,ea,sure of the relevant periodic process, normalised 
in such a way that n±{[k,k + 1] x T"^"^) = 1 for every k £ Z. Then there exist 
normalisation constants q± such that, 

hm {\f4A + fc) - q+^i+iA)\ + \n{A -k)- q-ii-{A)\) = . (5.4) 

A:— >oo 

(Here k is an integer.) Furthermore, this convergence is exponential, and uniform 
over the set A if we restrict its diameter. 

Remark 5.6. We used the shorthand notation A + k for {x + fc : a; G A}. 

Proof. We restrict ourselves to the estimate of ijl{A + fc), since the one on ijl{A — k) 
is similar. For fixed A; > 0, we introduce the sequence of stopping times given by 

4''^ = inf{t > : Xi(t) = fc} and then recursively a^n^ = ini{t > (fxn'^ : \Xi{t) - 
k\ — 1}, (fil^^i — inf{t > cri*^'' : Xi(t) — k}. This allows us to define an embedded 

Markov chain Z^^) on T'^''^ by setting zi'''^ = nX((/)^f^), where U{x,y) = y for 
{x,y) e M X T**-!. 

We similarly define an embedded Markov chain Z for the process X+. (By 
periodicity of X+, the choice of k is unimportant for the law of Z, so that we drop 
its dependence of k.) Denote by ir^^^ the invariant measure for Z^''^ and by tt the 
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invariant measure for Z. We then define cr-finite measures /Lt+ and \x^^^ on R x T''"-'^ 
through the identities 

^(fc) 

pP'\B)=l E,+fe,, / ' lB(X(s))ds7rW(cix) , (5.5) 

JT'i-i JO 

li+iB)= [ E^+fcei / ' lB{X+is)-k)dsTr{dx) . (5.6) 

Jt''-i Jo 

(Here and below we make a sHght abuse of notation and identify elements x G j^-i 
with the element (0,x) e K x T'*-^) It follows from [Has60, Thm 2.1] that ii^''^ 
is invariant for the process X and /i+ is invariant for X~^. Therefore, there exist 
constants > such that /x'*^-' = Cfe/x since the invariant measure for X is unique up 
to normalisation. Note that by translation invariance of X+, does not depend 
on k. 

Note that we can assume without any loss of generality that ^ C {x : xi > 0} 
(it suffices to shift it by a finite number of steps to the right in (5.4)). In this case, 
we can rewrite (5.5) as 

^W(A + fc)= / E^+feei / ' lA{X + {s)-k)dsTT''''\dx) . (5.7) 
JTd-i Jo 

This is because X{t) = X+{t) for t < a[''^ and, i{X{a[''^) < k, then /f,'/ 1a{X{s)- 

k)ds = whereas if X(crf ^) > k, then X{t) = X+{t) for t < This shows that 
the claim follows if we can show that ||7r — tt^*^^ ||tv — > as A; — > do and there exists 
a constant Coo such that ^ Cqo ■ 

Let lis first show that the latter is a consequence of the former. Setting Bk = 
[A;,fc + 1] xT"^-!, wchavecfe+i/cfe = ^W(Bfe+i)^('=+i)(Bfe+i). On the other hand a 

, (0) 

straightforward trial/error argument allows one to show that E^; j^^ '\.a{X~^{s)) ds 
is bounded uniformly over x € T''"^. It then follows immediately from (5.7) that 
there exists a constant C such that 

|/x('=)(Sfc+i) - M(Bo)| < C\\-K - ttW IItv , 
and similarly for \^^^'^'^\Bk^i) — ,u(i3o)|. It follows that provided that X]/j>o ll^r — 
tt^'^^IItv < 00, one does indeed have — >■ Coo- 

Denote now by P the transition probabilities for Z and by P'^^^ the transition 
probabilities for Z''^\ Then, wc can write P — QR, where R is the Markov kernel 
from T''-! to {-1,1} x T"^-^ given by R{x,A) = F^{X+{ai) G A) and Q is the 
Markov kernel from {-1, 1} x T''-^ to T''-^ given by Qix,A) = P^(X+((/)o) € A) 
for XiiO) = 0, (Ti = inf{i > : \Xi{t)\ = 1} and 0i = inf{t > ai : Xi_{t) = 0}. 
Since the diffusion X+ is elliptic, both Q and R are strong Feller and irreducible. 
It follows from the Doeblin-Doob-Khasminskii theorem [DPZ96, Proposition 4.1.1] 
that P{x, ■ ) and P{y, ■ ) are mutually equivalent for any x, y £ T''"^. Furthermore, 
it follows from the Meyer-Mokobodzki theorem [DM83, SeiOl, Hai07] that the map 
X P{x, ■ ) is continuous in the total variation topology. We conclude that the 
map (x, y) h- > |lP(x, • ) — P{y, ■ )|1tv reaches its maximum and that this is strictly 
less than 2, so that P satisfies Doeblin's condition. It follows that there exists a 
constant r] <1 such that P has the contraction property 

\\Pi/i - Pz/2||tv < vlWi - ^'2||tv , 
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for any two probability measures vi, vi on T ^ . Therefore, if we can find constants 
Efe such that 

sup \\P{x,-)-P^^\x,-)\\'v^<ek. (5.8) 

then we have 

IItt-tt^'^^IItv < ||P7r-P7r('=)||TV+||f7r('=)-P('=V('=)||TV < ?7||7r-7r('=)||TV+£fe , (5.9) 

so that IItt— TT^'^^ Ijxv < ^fe/(l ~ ??)• The problem thus boils down to obtaining (5.8) 
for an exponentially decaying sequence £fe. 

It follows from the same calculation as in Lemma 5.3 that the probability that X 
reaches the interface J?, before time when started on the hyperplane {x\ = A;} 
is bounded from above by ©(l/fc). This yields the "trivial" bound eu < ©(l/fc), 
which unfortunately isn't even summable. However, a more refined analysis allows 
to obtain Proposition 5.7 below, thus concluding the proof. □ 

Proposition 5.7. There exists a constant p G (0, 1) such that Sk < 0{p'^). 

Proof. The intuitive idea behind the proof of Proposition 5.7 is that if the process 

goes all the way back to the interface then, by the time it reaches again the plane 
{xi = k}, its hitting distribution depends only very little on its behaviour near the 
interface. In order to formalise this, let us introduce the Markov transition kernel 
(5+ from T''"^ to T''-! which is such that Q+{x,-) is the hitting distribution of 
the plane {1} x T'^-^ for the process X+ started at {0,x). Similarly, we denote by 
Q^''^{x, ■ ) the hitting distribution of the plane {k} x T"^"^ for the process X started 
at {i,x). 

For a fixed integer £ > r], our aim is to show that Q^''^{x, ■ ) gets very close to 
0+ ^ (a-', • )• Here, we denote by Q\. the A;th iteration of the Markov transition kernel 
Q+. With these notations at hand, define the quantities 

Q!fe= sup ||(9^'*^(a;,-)-(3+"^(a;,-)llTv , 

Note now that since, for fixed £, the probability that X reaches the interface 
before time (f)[''^ when started on the hyperplane {xi = k} is bounded from above 
by 0{l/k), we have 

ek< sup \\Q''-^^''{x,-)-Q+{x,-)\\ty (5.10) 

<5 sup \\Q^^'^{a:,-)-Ql-\x,-)\\TY<^ak, 

SO that it suffices to obtain an exponentially decaying bound on the a^'s. 

We now look for a recursion relation on the ctfc's which then yields the required 
bound. We have the identities Q^''= = 0'=-i'*=Q^ *-i and Q^'^ = Q+Q^'^-^. It 
follows from the triangle inequality that one has the bound 

WQ^'^S, - qI-'54tv < WiQ"-'-" - Q+)Q''''-'d4Tv (5.11) 

+ IIQ+(0'''-'4-Qt"'"''^-)llTv. 

At this stage, we note that by exactly the same reasoning as for P, the kernel 
satisfies Doeblin's condition. Therefore, there exists a constant fj <1 such that 

WQ+vi - (3+f2||TV < vlWi - ly^Wrv , 
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for any two probability measures vi, V2- This and the definition of Uk immediately 
implies that the second term in (5.11) is uniformly boimdcd by 770:^-1. On the 
other hand, it follows from (5.10) that the first term is bounded by ^oik, so that 

Oik < -j^Otk + mk-l , 
for some fixed constant C. The claim now follows at once. □ 

Finally, the last estimate that we need is the following. Denote by r the first 

hitting time of the interface and fix an arbitrary smooth positive function (p 
that is supported in the interval [1,2]. Set furthermore f^i^) = 'n~'^ip{n~^Xi) and 
(Pn{x) = n~^(p{—n~^xi). Then we have: 

Lemma 5.8. With the above notations, setting (p = ip{x) dx, we have 



[\tiX±it))dt-^ 

Jo -L'li 

uniformly for all x G {±ti} x T"^-! as n ^ 00. 

Proof. Again, we only consider the expression for X+ , the one for X_ follows in the 
same way. It follows from standard homogcnization results [BLP78, PS08] that the 
law of n^^X_|_(ri^t) converges weakly as n — > 00 to the law of Brownian motion with 
diffusion coefficient DJ*!- thus follows from [Bog07, Cor 8.4.2] that the law of 
n~^X+(n^t), where X+ is stopped at the first hitting time of ^ converges weakly 
as n — )■ 00 to the law of Brownian motion stopped when it reaches the hyperplane 

Denoting this limiting process by X"^, an explicit calculation allows to check 

that Ex JJ^ ip{X^{t))dt = ^ when xi = 1. Now, for any fixed T > 0, the map 

^11 

$T ■ X 1-^ (/?+(X(i)) dt is continuous, so that E^; /J" (p+(X+(f)) dt converges 
as n ^ 00 to J^^^ ip{X^{t)) dt. Letting T — )• 00 concludes the proof. □ 

We now have all the tools that we need to show that the exit probabilities from 
the interface converge to the desired limiting values. 

Proof of Proposition 5.1. Similarly to the proof of Proposition 5.5 we use a repre- 
sentation of the invariant measure fi in terms of an embedded Markov chain. This 
time, wc consider the stopping times ^q'^' = inf{i > : |Xi(t)| = 77} and then 
^1*=^ = inf{t > 4,i^^ : \Xi{t)\ = k}, 4>l^l^ = ini{t > ait^ : \Xi{t)\ = tj}. Denoting 
as similar to before by Tf^*^) the invariant measure of the embedded Markov chain 
Zn^ = X{^n^) (which is now a Markov chain on 9J^), we set 

il('^^B)= [ / ' lB(X(s))dsf«(dx) . (5.12) 

Jd.JK, Jo 



Again, the measures jl^''^ differ from /i purely through a scaling factor, so that there 
are constants Cfc such that il^''\B) = Ckl^iB) for every measurable set B. 

The idea now is to evaluate jl^^\ip^) in two different ways and to compare the 
resulting answers. First, we note from Proposition 5.5 that 
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On the other hand, combining Proposition 5.4 and Lemma 5.8 with the definition 
(5.12), we see that 

/x«(c/p±) = 4^+o(l) (5.13) 
^11 

as A; — > oo. Combining these two identities, we see that 

^ = ^ + 0(1), 

thus concluding the proof. □ 

6. Computation of the drift along the interface 

This section is devoted to the computation of the drift coefficients aj along the 
interface. Denote by r" the first hitting time of dJ^n by the process X. With this 
notation, recall that, by (4.4), we have the identity 

aj= lim -E^ / bj(Xs)ds, (6.1) 

n^oo n Jq 

provided that this limit exists and is independent (and uniform) over starting points 

Proposition 6.1. The expression on the right hand side in (6.1) converges to the 

expression given by (2.4). uniformly in x € J^ri- 

In order to show this, we will use the same construction as in the proof of 
Proposition 5.1. In particular, recall the definition (5.12) of the measures fl^''\ 
which are nothing but multiples of the invariant measure fi, as well as the sequence 

of stopping times ■* and ait^ . Denote furthermore by ni!'^ the invariant measure 
for the process on 9^ with transition probabilities P{x, A) given by 

P{x,A) '=>4X(0('=)) e A|r" > ^^t^) . (6.2) 

Our proof will proceed in two steps. First, we show that the limit (6.1) exists 
and is equal to the value (2.4) given in the interface, provided that we start the 

~ (k) 

process X in the stationary measure 7r„ and let /c — > oo. In the second step, wc 
then show by a coupling argument similar to the proof of Proposition 5.4 that the 
expression in (6.1) depends only weakly on the initial condition as n gets large, 
thus concluding the proof. 

Before we proceed with this programme, we perform the following preliminary 
calculation: 

Lemma 6.2. One has the normalisation 

l„t-*(«([-MlxT^-.)=2(|t + |^)SV, 

where the coefficients p± are as in (2.2). In particular, if n is normalised as in the 
introduction, then one has k~^jl^''^ w for large values of k. 

Proof. We know from Proposition 5.5 that iJ,{dx) — )• ii±{dx) at exponential rate as 
xi — > ±oo, so that on large scales /i behaves like a multiple of Lebcsgue measure 
on either side of the interface. Furthermore, we know from Proposition 5.1 that 
the corresponding normalisation constants satisfy the relation (2.2). Combining 
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this with the fact that jl'^^^ is just a multiple of fi, the result then follows from 
(5.13). □ 

Using this result, we obtain: 

Proposition 6.3. The limit 

1 r" 

a,- = lim lim — E.(fc) / bJXg)ds, 

[ {bj{x) + Cgj{x)) ,i{dx) , (6.3) 

where g is the function fixed in Section 3 and the constant /3 is as in Lemma 6.2. 

Remark 6.4. Note that if (p is any smooth compactly supported function, then 
the identity J £(j){x) fi{dx) = holds. As a consequence, the expression (6.3) is 
independent of the choice of the compensator g. 

(k) 

Proof. It follows from the definition of 7f„ and the strong Markov property of X 
that one has the identity 



exists and is equal to 



E.w r 6,-(X,)ds= 5^(P.(.)(^fUT"))'"E.(„ / ' ' h^{X,)ds 

" m>0 " " 

\,.,ff'-^^h,{Xg)ds 



Note now that it follows from Lemma 5.3 that 



(6.4) 



F{4>[^'' > t"") = k/n + 0{l/n) . (6.5) 

Since lim.n^ca 9jiX{T"'))/n = and furthermore, using the same argument as in 
(5.9), we have lim„^oo \\Ti'n'' — •^^'^-'Htv = for every > 0, so that 



lim -E (fc) / bj{Xs)ds= lim -E (k) [ bj{Xs) ds + gj{X{Tn 

1 r" ~ 

= lim — E.(fc) / hj{Xs)ds 

n^co n Jq 



)) 



1(1=) A n 

lim -rE^(fc) / bj{Xs)ds 

n^oo k Jq 

-E^(.) bj{Xs)ds 

bj{x) ii^''\dx) . (6.6) 



Here, we used (6.4) and (6.5) to go from the second to the third line and we used 
the definition of the /i'-'^-' to obtain the last identity. The claim now follows from 
Lemma 6.2. □ 

We can now complete the 
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Proof of Proposition 6.1. In view of Proposition 6.3, it remains to show that 
1 



hm 

n— foo n 



/ h{Xs)ds-¥.y I b{Xs)ds 
Jo Jo 



= 



uniformly over x,y G J^. Fix an arbitrary value of k > rj and consider again the 
transition probabilities P given by (6.2). Since they arise as exit probabilities for 
an elliptic diffusion, we can show again by the same argument as in the proof of 
Proposition 5.5 that P satisfies the Doeblin condition for some constant ri, namely 
\\Pi^i — Pi^2\\tv < (1 ~ 1)1^1 ~ uniformly over probability measures y\ and 

V2 on dJ^^. Note now that one has the identity 



r 6(X,)ds=^( W wSf^ <T^))W„, r ' b{X,)d8 (6.7) 

■'^ m>0 0<e<m •'^ 

m>0 



where we denote by (resp. E^) the probability (resp. expectation) for the process 
X started at P'"(a;, •)• 

Note now that we have the identity 

P.(^^' < r") = P. (^f ) < T") + Pp.(,,. ) {^^^l, < T") . 

Also, by choosing k sufficiently large (but independent of n), we can ensure that 
there exist constants c, C > such that 

l--<P.(^l'^<r")<l--, 
n n 

imiformly for x G .X; and for n sufficiently large. It also follows from the contraction 
properties of P that 

< r") - PL(0i'^ < r")| < 2(1 - nr , 

uniformly over x,y G ^j^. 

Combining these bounds, we obtain for every l<m l\n the estimate 

K ( 

|Px(0L'^ < r") - P,(0« < r")| < — + 2(1 - . 
In particular, there exists a constant K, such that wc have the uniform bound 

valid for every m > and every n sufficiently large. Summing over m, it follows 
that 

^ |P,(0W < r") - P,(0W < r")| < i^v^ , 

m>0 

for a possibly different constant K. 

On the other hand, it is possible to check that there exists a constant C (de- 
pending on k) such that 

if")/ 



'0 



[<P\ At 

/ h{X,) ds 

Jo 
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uniformly over a; G J^, so that 

E: 



b{Xs) ds-W^ [ ' b{Xs) ds < 2C(1 - r])'' 
Jo 

Inserting these bounds into 6.7, we obtain 

/ b{Xs) ds-Ey f b{Xs) ds < 2C V (1 - r/)" + C^/E , 
•^0 -^0 m>0 

so that the requested bound follows at once. 



□ 



6.1. Bound on the second moment. In order to conclude the verification of the 

assumptions of Theorem 2.4, it remains to show that the second bound holds in 
(4.4). For the non-rcscalcd process, we can reformulate this as 

Proposition 6.5. For every fj > 0, there exists a constant C > such that the 
bound 

E,||y(T")-2/f <Cn^ 
holds for every n > 1 and every initial condition y G J^fj. 

Proof. It follows from (3.1) that 

2 



Ej, F(t") -2/ < 2Ej 



b{Xs)ds 



a{X,)dWis) 



It follows from Ito's isometry that the second term is bounded by CEt", which in 
turn is bounded by 0{n^) by a calculation virtually identical to that of Lemma 4.5. 

It remains to bound the first term, which we will do with the help of a decompo- 
sition similar to that used in the proof of Proposition 5.1. For two constants c > 
and a > to be determined, we set cpo = 0, o-„ = inf{t > '■ \Xi{t)\ = c + a}, 
and ^„ = mi{t > an-i ■ \Xi{t)\ = c}. Define furthermore 

N = inf{fc > : (7fe > r"} . 

Since b is supported in a bounded strip around Jf), we can make c sufiiciently large 
so that the first term in (6.8) is bounded by some multiple of 



E, 



N 

•y i^i^k 



fe=0 



< 



N 

E„Af3Ey^(afe-</.fe)4 



fe=o 



< 



EyN^J^Mi'^k - <kkY I N > k)fy{N > k) . 



fe=0 



Note now that since ak is the exit time from a compact region for an elliptic 
diffusion, there exists a constant C such that Ej,((c7fe — (/>fc)* | A'^ > fc) < C, uniformly 
in y. Furthermore, it follows from Lemma 5.3 that if a is sufficiently large, then 

fy{N >!)<!--, 
n 

for some constant c > 0, uniformly in y. The strong Markov property then imme- 
diately implies that W'y^N > k) < (1 — ^) , so that N is stochastically bounded by 
a Poisson random variable with parameter 0{n) and the claim follows. □ 
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7. WeLL-POSEDNESS of the martingale problem and CHARACTERIZATION 



The aim of this section is to show tliat the martingale problem associated to 
the operator C as defined in Theorem 2.4 is unique and to characterise the corre- 
sponding (strong) Markov process. Our main tool is the following general result by 
Ethier and Kurtz [EK86, Theorem 4.1]: 

Theorem 7.1. Let E he a separable m,etric space, and let A: ^{A) — >■ Bb{E) be 
linear and dissipative. Suppose there exists A > such that 



and such that C is separating. Let ji G ^{E) and suppose X is a solution of 
the martingale problem for {A, /x) . Then X is a Markov process corresponding to 
the semigroup on C generated by the closure of A, and uniqueness holds for the 

martingale problem, for [A,fj,). 

See also [BP87] for a more general result on the well-posedness of a martingale 
problem with discontinuous coefficients. This allows us go finally give the 

Proof of Theorem 2.4- Since we already know from the results in the previous two 
sections that limit points of X^ solve the martingale problem associated to £, it 
suffices to show that this martingale problem is well-posed and that its solutions 
are of the form (2.6). 

For this, we somehow take the reverse approach: first, we construct a solution 
to (2.6) and we show that this is a Markov process solving the martingale problem 
associated to £. We then show that this Markov process generates a strongly 
continuous semigroup on Co(K'^), whose generator is the closure of £ in Co- Since Cq 
is separating and since generators of strongly continuous semigroups are dissipative 
and satisfy (7.1) by the Hille-Yosida theorem, the claim then follows. 

In order to construct a solution to (2.6), let M± be matrices satisfying M±M^ = 
and such that 



for some vectors v± G M**"^ and some (d — 1) x (cZ — 1) matrices M±. (This is always 
possible by the QR decomposition.) We then first construct a Wiener process Wi 
and a process Xi such that 



where L is the symmetric local time of Xi at the origin. This can be achieved for 
example by setting Xi = g{Z), where 



OF THE LIMITING PROCESS 



C =^^{X-A) = @{A) , 



(7.1) 






Z is a skew-Brownian motion with parameter 
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and W is the martingale part of Z. Given sueli a pair {Xi,W), we then let W 
be an independent d — 1- dimensional Wiener process and we define pathwise the 
M''~ ^-valued process X by 

X{t) = f\lx,<oM- + lx,>oM+) dW{t) + [\lx,<oV- + lx,>oV+) dW{t) 
Jo Jo 

+ a dL{t) , 
Jo 

where aj = ctj+i. Since we know that skcw-Brownian motion enjoys the Markov 
property, it follows immediately that Xi is Markov, so that X = {Xi,X) is also a 
Markov process. Applying the symmetric Ito-Tanaka formula to f{X) it is further- 
more a straightforward exercise to check that X does indeed solve the Martingale 
problem for C. 

The corresponding Markov semigroup {Vt}t>o maps Co(M'') into itself as a con- 
sequence of the Feller property of skcw-Brownian motion [Lej06]. Furthermore, as 
a consequence of the uniform stochastic continuity of X, it is strongly continuous, 
so that its generator must be an extension of C. Since the range of jC contains 
Cq°(]R''), which is a dense subspace of Co{M.'^), the claim follows. □ 
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